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in the calculation. On the other hand, we show that Beyond the Standard Model physics 
with a characteristic scale comparable to the scale that governs the homogeneous decay 
rate in the SM, can in principle lead to an enhanced decay rate. 

ArXiv ePrint: 1509.05405 



Contents 


1 Introduction 1 

2 Derivation of the Heterogeneous Decay Rate 3 

3 Thin-Wall Examples 8 

3.1 Scalars 9 

3.2 Fermions 10 

4 Application to the Higgs Vacuum 11 

5 Discussion and Conclusions 14 


1 Introduction 

The discovery of the Higgs boson combined with the absence of any observation of physics 
beyond the Standard Model (SM) has an intriguing implication for the ultimate fate of the 
universe. With the measured value of the mass of the Higgs boson as the final SM input, 
m,H = 125.09 ± 0.21(stat.) ± 0.11(syst.) GeV [1], state of the art calculations suggest that 
the electroweak (EW) vacuum is unstable [2-4] (see also references therein). However, the 
lifetime of the universe is computed to be many, many orders of magnitude larger than 
the current age of the universe; the EW vacuum is said to be metastable. Due to this 
metastability, there is no need to invoke beyond the Standard Model (BSM) physics to 
explain the observed age of the universe. See [5] for a recent review. 

Much of the intuition developed for calculating the decay rate of an unstable phase of a 
quantum held theory (QFT) is based on theory developed for systems that can be studied 
in the laboratory. In fact, the framework developed by Coleman and collaborators [6— 
8] for calculating decay rates in QFT, which is what is used to determine whether the 
Higgs vacuum is metastable or unstable, is the relativistic, four-dimensional analog of the 
previously developed methods used in statistical mechanics. 

In particular, what was calculated when it was suggested that the universe is metastable 
is the homogeneous decay rate. However, in most commonly studied systems it the hetero¬ 
geneous decay rate (rather than homogeneous) that dominates the catalysis of the phase 
transition. In a typical laboratory system, a phase transition is seeded by the presence of 
an impurity either in the bulk of the unstable phase, or by the boundary of the system. 
Since the intuition for the Higgs vacuum decay rate has been developed from these systems, 
and since the universe is not just a constant electroweak vacuum, it seems worthwhile to 
study the heterogeneous decay rate of the Higgs vacuum. 
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Classical nucleation theory, including the heterogeneous case, seems to have first been 
introduced by Volmer for liquid nuclei [9, 10]. The extension to the case of crystal nucle¬ 
ation on surfaces was made by Turnbull [11-13] (and their respective collaborators). The 
nucleation theories that most closely resemble Coleman’s work were subsequently developed 
by Langer [14, 15]. 

There is a relatively small amount of literature on induced vacuum decay in QFT. 
Refs. [19, 20] looked at the induced decay rate due to single particles. An analysis based 
on these methods, Ref. [21], was done to see if cosmic rays could catalyze vacuum decay. 
Classical catalysis, rather than quantum tunneling, due to Hawking radiation has been 
investigated in Refs. [22-24], Black holes as the seeds of bubble nucleation were discussed 
in Refs. [25-31]. 

Enhanced vacuum decay rates have been investigated outside the context of an impu¬ 
rity induced decay as well. A general analysis of tunneling in theories with multiple scalar 
fields was made by Refs. [32, 33] using the thin-wall approximation. The decay rate of the 
Higgs vacuum at finite temperature was recently updated in Ref. [34], While Refs. [35-41] 
investigated the stability of the Higgs vacuum in the early universe with a particular focus 
on inflation. 

When studying vacuum decay outside of the SM, it is natural to expect that the 
lifetime of the electroweak vacuum would in general be different from the prediction of the 
SM. For example, Refs. [42-44] colorredconsidered the SM supplemented by Planck-scale 
suppressed higher-dimensional operators such that the vacuum structure of the theory is 
different from that of the SM. In was shown in [42-44] that the lifetime of the EW vacuum 
in these theories can indeed be much shorter than it is in the SM. For additional discussion 
of this scenario, see [5, 41]. 

The rest of the paper is organized as follows. In Section 2, we derive the decay rate 
of an unstable phase of a quantum field theory in the presence of an impurity in the thin- 
wall approximation. This derivation is based on the how the impurity changes the (flat 
spacetime) geometry relative to case of pure false vacuum. Then, in Sec. 3, two examples 
are given that show how to estimate some of the additional parameters that enter into this 
heterogeneous decay rate. After that, Sec. 4 applies the formalism developed in previous 
two sections to the EW vacuum of the SM. The goal is to determine if the conclusion 
drawn in the homogeneous analysis about BSM physics not being necessary to explain the 
observed age of the universe is still valid in the heterogeneous case. We find that baryonic 
matter (stars), which acts as an impurity in the electroweak (metastable) vacuum, leads 
to a heterogeneous decay rate that is suppressed with respect to the homogeneous case. 
That is to say, the conclusions drawn from the homogeneous case are not modified by the 
inclusion of baryonic matter in the calculation. Additionally, we confirm that BSM physics 
with a characteristic scale comparable to the scale that governs the homogeneous decay 
rate in the SM (or very dense physics) can lead to an enhanced decay rate. 

1 This is likely an incomplete list of references as the subject is rather old. However, as a consolation, 
see [16-18] for some reviews of the experimental confirmation of classical nucleation theory and data on the 
decay rates of superheated liquids. 
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2 Derivation of the Heterogeneous Decay Rate 


In this section, we derive the decay rate of a metastable phase in the presence of an impurity 
in the thin-wall approximation. The decay rate can generically be written as, 


l= A e~ B , ( 2 . 1 ) 

where S' is a (in general dimensionful) symmetry factor to be discussed in what follows. In 
QFT, the coefficient B is equal to the change in the Euclidean action of the system due to 
the appearance of a bubble of stable phase within the bulk metastable phase, B = S^. In 
the thin-wall approximation, this change in the action is given by, 

Se, i = - Y r i e i3 + 2vr 2 rf5i 3 , (2.2) 


where we are using a notation similar to that of Coleman [ 6 ] (the subscripts are for later 
convenience) with 633 and S '13 being the difference in the energy density and the surface 
tension between the true and false vacua respectively, ri is the radius of the bubble. 
Extremizing the action with respect to the radius yields the critical bubble size needed for 
nucleation to occur, 


r* = 


3^13 

£13 


Substituting r\ back into Se yields the critical action, 


(2.3) 


5* 


EA 


27t r 2 Sf 3 

2e? 3 


(2.4) 


In this work, we will only be interested in order of magnitude estimates for the pre¬ 
exponential factor, A, but it is an important factor to keep track of nonetheless. To be 
more precise regarding Eq. (2.1), the quantity of interest for a homogeneous QFT is the 
decay rate per unit volume, F/V 3 . Recall that this is so because the coefficient B is 
invariant under (Euclidean) spatial translations in this case, leading to the probability for 
tunneling being proportional to the volume of spacetime. Then, by dimensional analysis, 
the pre-exponential factor of T/V 3 must be proportional to l/(r£) 4 . 

In the case of heterogeneous nucleation due to the presence of an impurity in the bulk, 
the bubble is no longer a full sphere. Instead, the generalized bubble has a lenticular shape 
formed by joining two spherical segments as shown in Fig. 1; see Fig. 2 for a couple of 
special cases. It is this decreased volume of the bubble in the heterogeneous case that 
leads to the smaller action, and thus in principle a faster decay rate. We assume in this 
analysis that the characteristic size of the impurity is much larger than that of the bubble 
such that the boundary between regions 1 and 2 can be treated as a flat plane. It would 
be interesting to relax this assumption to study the cases where the impurity is a string 
or monopole of comparable size. Fig. 1 is a two-dimensional projection of the bubble, and 
is symmetric about the horizontal (blue-red) line in the additional (Euclidean) dimensions 
that aren’t shown. 
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Figure 1: Two-dimensional projection of a bubble of the stable phase (region 3) that has 
formed on the boundary of a bulk metastable phase (region 1) and an impurity in the 
bulk (region 2). The characteristic size of the impurity is assumed to be much larger than 
that of the bubble such that the boundary between regions 1 and 2 can be treated as a 
flat plane. This picture is symmetric about the horizontal (blue-red) line in the additional 
dimensions that aren’t shown. 



Figure 2: Special cases of the situation presented in Fig. 1. In Fig. 2a, a full sphere forms 
within the impurity. While in Fig. 2b, a “half-bubble” forms on the interface of regions 1 
and 2, but only on the side of the impurity. Note that the distance between the bubble 
and the interface between phases 1 and 2 is not necessarily to scale in Fig. 2a. See the text 
for more details on these special cases. 
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The Euclidean action in the heterogeneous case can be written as follows, 

S E ,e = S E , i f(0i) + S E , 2 f (0 2 ) - ^vrr 3 (S 12 - S&) , (2.5) 

7T 2 

S E ,i = ——r(e,; 3 + 2n 2 rfSi3, 

5^2 = S’l3 COS 91 + 523 COS 9 2 , 

where the subscript 6 labels the heterogeneous case. In the preceding equations, e^- and Sij 
are the energy density difference and the surface tension between phases i and j respectively. 
As in Fig. 1, regions 1, 2, and 3 correspond to the metastable phase, the impurity, and the 
stable phase respectively. S) is action associated with a bubble that forms completely in 
region i, see Eq. (2.2). The form of the first two terms on the right-hand side of Eq. (2.5) 
come from cutting off the angular integration of the spheres at 0y 2 . / (9) is the fractional 
volume of a 4-sphere capped at angle 9 rather than it (i.e. / (7r) = 1, / (n/2) = 1 / 2 , etc.). 
The relevant geometry results can be found in e.g. [45]. The last term on the right-hand 
side of Eq. (2.5) comes from change in the action caused by the appearance of the stable 
phase, which eliminates some of the hyperplanar boundary between the metastable phase 
and the impurity . 2 Lastly, 5 32 is the value 5 32 takes when a bubble of critical size forms on 
the interface between the metastable phase and the impurity. As can be seen from Fig. 1, 
only three of five length/angle parameters in Eq. (2.5) are independent. One way to write 
the constraints is, 


r = r\ sin 9\ = r 2 sin 9 2 . 

From this, it immediately follows that, 

(2.6) 

h\,2 = n, 2 (1 - COS0i i2 ) . 

(2.7) 


As was done in the homogeneous case, we now look for the critical points of the 
heterogeneous action, using r, h±, and h 2 as our free parameters. However, in what follows 
we will switch between various combinations of r, hi, ri, and 9i to write the equations in 
their simplest form. We find 12 unique critical points, five of which can be eliminated by 

2 There are an infinite number of ways to return to the homogeneous action from the heterogeneous case. 

A particularly simple way is to send —> 7r; 62 , r —> 0. 


f(0) = - 

7r 


9 — sin 9 cos 9 1 + 
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requiring r, h\ . and h 2 to be real and non-negative. The first derivatives of Se,o are, 


dSE,e 

dr 

dSE,e 

dhi 


Se3,i 9 (@i) + Se3,2 9 (# 2 ) — 4vrr 2 (5 i 2 — S * 2 ) , 
irrf (~ri€i 3 + 35 i3 ) h (6i ), 


Se3,i = -^irrfe i3 + 4vr r 2 5* 3 , 

g (6) = - cot ^(6 9 — 3 sin ( 9 ) — 3 sin (29) + sin (39)), 
9 cos (9) — ^ sin (0) — sin (39) J . 



( 2 . 8 ) 


The simplest critical points to analyze are when a full sphere (or no sphere) forms 
completely within the bulk of the metastable phase and/or the impurity, r* —>■ 0 ,/?,*—>■ 
{0, 65 * 3 / 6 * 3 }. For these cases, the critical action takes simple forms, 


5* 


Efi 


— (o> 


s* 


E, 21 


$E ,1 + ‘Sz?,2}- 


(2.9) 


The first solution is the case when no bubble forms, which we aren’t interested in. While 
the second solution is just the homogeneous case. The third solution is a bubble that forms 
completely within the impurity, as in Fig. 2a. For some range of parameters, this solution 
will be energetically favorable with respect to the homogeneous case. The same cannot be 
said for the last of these solutions, which is always energetically disfavored compared to 
the homogeneous case. 

In addition to being energetically favorable with respect to the homogeneous case, for 
a heterogeneous solution to be of interest, the bubble must grow. The second derivatives 
of Se, 9 , evaluated at the critical values for r*, are, 


d 2 S E ,> 


dK 2 


d 2 SE,e 


r ;=0 


dr 2 


= 0 , 


( 2 . 10 ) 


n= 0 , r2=0 


d 2 SE,e 


dh ? 


3S i3 

r -LT 

e i3 


d 2 SE,e 


dr 2 


ri= 4w’ rj=0 


— cos (9i) CSC 2 ( \ h ( 9i ), 

- 37 r— cot 2 ( ^ (6 9i - 8 sin (<9*) + sin (20*)), 
V 4 / 


d 2 SE,e 

dr 2 


3Sl3 *. 3^23 

e 13 £ 23 


d 2 SE,e 

dr 2 


ri = 3513 , r2=0 
e 13 


+ 


d 2 SE,e 

dr 2 


r2= ^ ,ri=0 


For the third solution of Eq. (2.9), the second partial derivative test is ambiguous because 
one of the eigenvalues is zero, whereas the last solution is found to be a local maximum. 
However, a numerical study shows that the second and third solutions in Eq. (2.9) are 
saddle points of Se,o (with respect to r, hi, and h 2 ). Thus, there is the possibility of 
bubble growth after formation. All the cases that are discussed in what follows are also 
saddle points. 
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The dynamics of bubble growth in Lorentzian space are easy to picture for these 
solutions. The bubble grows similarly to how a bubble would grow in the homogeneous 
case, as described by Coleman [ 6 ]. The only difference being that for a bubble forming 
completely within an impurity, the bubble first converts the impurity to true vacuum, 
before continuing to expand (at the speed of light) out into bulk false vacuum, which is 
also then converted to true vacuum. 

The second most simple case is when a bubble forms at the interface between the 
metastable phase and the impurity, but only on one side of the interface as in Fig. 2b. We 
dub this scenario the “half-bubble.” An example of such a scenario is briefly discussed in 
Sec. 3.1. First consider the case where h -2 = 0. The case where hi = 0 is trivial to obtain 
from the solution with I 12 = 0. The critical sizes in this scenario are, 

3 (S 13 — S 12 + S 23 ) ^\Z^i3 — (*^ 12 — ^ 23 ) 

£13 ’ £13 

From this, we see that S 12 — S 23 = S 13 cos 9\. The critical action in this case is, 




$E,e 


= s* Eil net) 


( 2 . 12 ) 


with cos 9* = (S 12 — S 23 ) /S 13 . To obtain the analogous solution for when a bubble forms 
on the interface within the impurity, but not in the metastable phase, simply switch the 
l’s and 2 ’s in Eqs. ( 2 . 11 ) and ( 2 . 12 ), i.e. Sg g = 2 / (# 2 ) with cos 9% = (S 12 — £ 13 ) /S 2 3 . 

The remaining solution is the most general scenario (in terms of parameters), where a 
bubbles forms at the interface and on both sides of the interface, see Fig. 1. In this case, 
it is easier to express the conditions for criticality in terms of r* and 9\, rather than r and 
K, 



(2.13) 


for i = 1,2 and j = 2,1. Note that there is a similar critical point, with a plus sign rather 
than minus sign in its solution for cos 9*. However, this leads to an imaginary value for r; 
cos 9i = 5'i2£i3 + 70 r 2 < 0. Just as in the half-bubble case, S 12 = S^ 2 , leading to a 
critical action of the following form, 


$E,e — 1 f {0 1) + Se ,2 f ($2) • ( 2 - 14 ) 

It will prove useful to examine the limit ej 3 —> 0, in which the angles take the following 
forms, 

cos 9i = 12 ——, cos 9j = 1. (2.15) 

*->i3 
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These angles are the same as in the half-bubble case, and in fact this limit reduces to 
the half-bubble case. We now a have way to determine which side of the interface the 
half-bubble will form on, and that is the side i, on which e *3 7 ^ 0 . 

Lastly, we discuss the pre-exponential factor in the heterogeneous case, Ag. This factor 
can be thought of as an attempt frequency that takes into account the numbers of sites 
available for a heterogeneous decay to occur [46]. Assuming all the impurities are the same, 
in analogy with the homogeneous case, quantity of interest in the heterogeneous case should 
the decay rate per number of impurities per unit surface area, F/(Ni mp S 2 )- If this is so, 
then Ag should be proportional to l/(r *) 3 by dimensional analysis. Because there is a 
smaller amount of symmetry in the heterogeneous case, Ag is typically smaller than A, 
meaning that Se,o cannot be arbitrarily smaller than Se and still increase the decay rate . 3 
There is some minimum decrease in the action needed to overcome the suppression of the 
pre-exponential factor to obtain an enhanced decay rate. 


3 Thin-Wall Examples 


In the previous section, we showed that the action associated with the decay rate of a 
metastable phase in the presence of an impurity has several additional parameters with 
respect to the homogeneous case. The question then becomes, how does one calculate those 
parameters from a microscopic theory. As it is the simplest case, we will mostly focus on 
the scenario where the bubble forms completely within the bulk of the impurity, Fig. 2a. In 
particular, we give two examples in the thin-wall approximation that show how to calculate 
S 23 and 623 in these examples. In both examples, we give an expression for Se,q, but save 
the analysis of the heterogeneous decay rate for Sec. 4. We hope that these examples will 
serve a guide to those who wish the compute the heterogeneous decay rate other models. 

Before we get into the examples, we give a quick review of the thin-wall model of [ 6 ]. 
The Euclidean Lagrangian is, 

£e = 2 (<9</>i) 2 + (0i “ a i) + 2 ^ - ai ^ ’ ( 3 -l) 

where the subscripts are for later convenience. The solution, 4>i(p), is invariant under four¬ 
dimensional Euclidean rotations, with p = \Jt 2 E + |x | 2 being the 4 -d radial coordinate. The 
thin-wall approximation is X\af 613 , such that the energy difference between the two 
vacua is small, and that the Lagrangian has a Z 2 symmetry to an excellent approximation. 
In the thin-wall approximation, cf> 1 can be expressed as, 


01 (P) 


-a 1 

< a± tanh (p - ri)) 
o\ 


p <C r\ 

P ~ ri , 

p > n 


(3.2) 


with pi = y/X[ai and r\ being the radius of the bubble (in the homogeneous case). From 
Eq. (3.2), we see that the field configuration deep within the bubble is the true vacuum, 

3 More precisely, the pre-exponential factor in the probability (rather than decay rate) is typically smaller 
in the heterogeneous case. 
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and that far outside the bubble is the false vacuum. While near the wall, Eq. (3.2) shows 
that the field configuration is a soliton, which smoothly interpolates between the two vacua. 
This solution for 0i yields a Euclidean action of the form of Eq. (2.2) up to corrections of 
order O (^ 0 2 r]“ 2 ), which are small in the thin-wall approximation. With the solution for 
01 . the surface tension of the bubble can be calculated in the usual way, 

r°° 9 9 

<Si3 = J dpCs = -\Ai«i = (3-3) 

There are a couple of changes in the heterogeneous case. One difference is that, in 
general, the impurity breaks the 0(4) symmetry of the homogeneous solution. Non-0(4) 
symmetric tunneling solutions have been studied by Refs. [25-28] in the context of black 
holes as the seeds of the bubble nucleation. In addition, it is well known that the tun¬ 
neling solution at finite temperature possess an 0(3) symmetry, rather than the full 0(4) 
symmetry [34, 47]. 

Another notable feature of the heterogeneous case is that in general there are additional 
field configurations for 0i, beyond those listed in Eq. (3.2). These regimes are due to the 
impurity, and can be found by solving the equation of motion for 0l including the 0i field’s 
interactions with the impurity. Regrettably, we do not know of any such analytic solutions 
for the examples we discuss in what follows. However, some progress has recently been 
made in this direction, at least in the case of radiatively generated potentials [48, 49], To 
proceed, we do not consider the full parameter space of these examples, but rather we 
make choices for some subset of parameters that allows to use the homogeneous solution 
for 0i in our calculation of the heterogeneous decay rate. In any case, S 23 is given by an 
expression analogous to Eq. (3.3) except that the Euclidean Lagrangian now also includes 
0i’s interactions with the impurity. Typically, 623 can simply be read off of the Lagrangian 
that includes 0 i’s interactions with the impurities, as is done for £13 in the homogeneous 
case. 

3.1 Scalars 

As a first example, consider adding a second scalar held to the thin-wall model of [ 6 ]. The 
Lagrangian is, 

£ = \ (<90i) 2 + ^ (<90 2 ) 2 - U (0i, 02 ) , (3.4) 

U = ^ (0? - a?) 2 + y (02 - a!) 2 + y (</>i - a?) a 2 (0 2 - a 2 ) + (0i - 01). 

Here 0i is the bounce solution that governs the tunneling rate, and 0 2 is the additional 
scalar held. We will assume the held configuration of 0 2 is that of a kink, say, in the z- 
direction, centered about z 0 , which breaks the 0(4) symmetry of the homogeneous solution 
and provides a preferred direction for tunneling. We will also assume that the surface 
tension of the 0 2 kink is large enough that the backreaction of 0 i on 0 2 is negligible, and 
that it is narrow enough such that the boundary between regions 1 and 2 can be treated 
as a hat plane. Both of these assumptions can be accommodated by taking /r 2 = \f\ 2 a 2 to 
be much larger than any other energy scale in the problem. 
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For z > zo, the A 3 interaction has no effect as (j) 2 ~ 02 in this region except for the 
transition near the kink of fa. However, this transition region should be small if the kink 
appears as a plane to the bubble, which we have assumed. Therefore, the surface tension in 
this direction, S13, is the same as in the homogeneous case. For z < zo, the effect of the A 3 
interaction is to shift the mass term of the first scalar field, [i\ —> fi\ + A 30 ^, where again, 
we have neglected the width of the kink of fa. The surface tension of the wall between the 
bubble of stable phase and the impurity is then given by S 13 , but with the aforementioned 
shifted mass, 

5 23 «^(^ + A 3 a |) 3/2 = 5 13 ( 1 + ^|) . (3.5) 

The difference between the energy density of the two vacua now depends on the shifted 
mass as well, 

£23 = e i3\/l + t"~—§• (3.6) 

V Aiaf 

From this, we see that for A 3 > 0, there is no possibility of an enhanced tunneling 
rate for a bubble forming completely within the impurity, as the action increases in this 
scenario. However, when A 3 < 0 and |A 3 | a| < Aiaf, the heterogeneous action for a bubble 
to form within the impurity is given by, 

s « = Se (' ■ • (3 - 7) 

which does lead to the possibility of an enhanced decay rate contingent upon the change 
in the pre-exponential factor in the heterogeneous case relative to the homogeneous case. 

On the other hand, for A 3 < 0 and | A 3 1 a| > \\a \, there is a unique vacuum at fa = 0 
for z < zq (ignoring the tiny width of the fa kink) with an energy density of Aiaf /8 up 
to corrections of order € 13 . While for z > zo, the two vacua of the homogeneous case are 
still present at fa = ±a 1. Comparing energy densities, we see that £23 = O(Aiaf) ei 3 . 
Based on this comparison and the discussion of Sec. 2 , the picture for this scenario is 
given by Fig. 2b. However, this comparison also calls into question whether the thin-wall 
approximation is still valid, as 623 ~ S 13 . For this reason, we do not consider this range of 
parameters further in this work. 

3.2 Fermions 

Consider adding a fermion to the standard thin-wall example of [ 6 ]. The fermion has a 
Yukawa interaction with the scalar, 


U D yi>ipfa. (3.8) 

Assume that in some regions of space, there is a non-zero fermion number density, (|x|) = 
('ijjfa 7 ^ 0. These fermions might condense into composite objects, like in QCD, or there 
may simply be a non-zero number of elementary fermions. In either case, we will assume 
that the scalar interacts with this object through the Yukawa term. 
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The difference between the homogeneous and heterogeneous cases is the addition of the 
Yukawa interaction, so 4 >l will be required to tunnel through the finite density of fermions in 
some regions of space. The surface tension between the metastable phase and the impurity 
should be given by, 

S23 = Si3 + yN^[ dp<f>i(p), (3.9) 

Jo 

if the following two assumptions are satisfied: (1) L ?'i i2 , in order to satisfy the assump¬ 
tion that the boundary between the metastable phase and the impurity appears as a plane 
to the bubble, and (2) the fermion density is constant, Np (|x|) = N^, in the regions of 
space where it is non-zero. Experience from the real world suggests that L > ' ; the 

majority of the nuclei in the universe are light, and they are all at least as big as a proton, 
which has a size ~ -^qcd ■ Taking L >, A/, 7 , the planar assumption suggests that Np 
is small enough such that the homogeneous solution for is not significantly perturbed, 
which yields 

S 23 - Si 3 « In (cosh V /3 )) ' (3-10) 

Combing the planar assumption (I > n) with the thin-wall approximation (Aiaf ei 3 ), 
we find y/X\aiL ^$> Xiaf/eis 3> 1. Thus, Eq. (3.10) can be rewritten as, 

S 23 - Si 3 « yN 2 J \i - 0{Np). (3.11) 

The difference in the energy density between the false and true vacua in the heteroge¬ 
neous case can be read off of the Lagrangian, 


£23 — £i 3 + 2yA^ai. 


(3.12) 


For a bubble that forms completely within the impurity, the critical Euclidean action in 
the heterogeneous case is given by, 


Se,9 


Th t 2 ( Si:i + y N l /3(i i) 
2 (ei 3 + 2yNpai) 3 


Se 



4 yNfW 

Si 3 



(3.13) 


For y < 0, the preceding equation gives a suppressed action, and thus the possibility of 
an enhanced decay rate. However, since JVh has been assumed to be small, its effect on 
the action will be small in this example regardless of the sign of y. Note that to decrease 
the action with these assumptions, y must be negative whether —a or a (ei 3 —> — ei 3 in 
Eq. (3.1)) is chosen as the true vacuum in the homogeneous case. 


4 Application to the Higgs Vacuum 

In this section, we apply the formalism developed in the previous two sections to the case 
of the Higgs vacuum with baryonic matter as impurities, as well as investigate generic BSM 
physics. 
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We start by briefly reviewing the homogeneous vacuum decay rate and probability. 
See Refs. [5, 50] for more details. In the case of the Higgs vacuum, the decay rate per unit 
volume in the SM is given by 


L ~ J_ e ~S{ Ab) 

v 3 R 4 


(4.1) 


Since quantum corrections enter the potential logarithmically, one must go to large field 
values to determine if there is an instability, h v. With this approximation, the bounce 
action in Eq. (4.1) can be obtained using the tree level potential, 


S(A b ) 


8vr 2 

3 |A (Ab)| 


(4.2) 


This action is classically scale invariant, and the size of the bounce, R , is arbitrary. Quan¬ 
tum corrections break scale invariance, contributing to action as AS ~ log (RAb)- Choos¬ 
ing A B R 1 minimizes these corrections, and resolves the implicit ambiguity in choosing 
which MS scale to evaluate A at. As always, one is interested in the bounce that maximizes 
the decay rate. With this choice of scale, the maximization condition is in practice given 
by (3\(Ab) = 0, where A b ~ 2.0 x 10 1 ' GeV in the SM [3]. Eq. (4.1) can be integrated over 
the past light cone of the universe to give the probability that the Higgs vacuum decayed 
at some point in the past light cone of the universe, 

p ° * Gt) ^ 3<Ab) ’ (4 - 3) 

where Hq is the Hubble constant. Numerically, we find Pq ~ 10 741 , compare against the 
left panel of Fig. 7 of Ref. [3]. 

After clarifying the assumptions that go into the calculation, we proceed to investigate 
the case of baryonic matter (i.e. stars) in the universe acting as impurities to seed the 
decay of the Higgs vacuum. 

First of all, it should be noted that the examples of the previous section relied on the 
thin-wall approximation, which is not valid for the SM. One way around this problem, 
which was the approach of [31], is to include additional, Planck-scale suppressed operators 
(for example: (WH) 3 /M^) to modify the Higgs potential such that it does satisfy the 
conditions for the thin-wall approximation to be used. We do not include any such oper¬ 
ators in our calculation, as we believe our results will not qualitatively be affected by this 
approximation, and hope that it will be self-evident as to why we think so in what follows. 
Furthermore, if such operators were included, it may be difficult to separate the effects of 
the impurity from the effects of the operators. 

In addition, it was shown in Sec. 3.2 that (at least for a bubble forming completely 
within the impurity, and a\ <C S 13 ) the Yukawa coupling of the fermions to the scalar 
field must be negative to actually cause a decrease in the action. This is also not the case 
in the SM, and there are at least two ways out of it. The first is to simply choose a more 
complicated critical point of the heterogeneous action. With a larger set of parameters, it 
will allow for a decreased action with a positive Yukawa coupling. A second choice is to add 
an operator (or two) to the theory of the form, (W H)(qid\H) or (H^H)(qiU\H), where 
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the subscripts are generational indices. These operators can make the effective Yukawa 
couplings of the up- and down-quarks negative, while keeping their masses positive. The 
first generation Yukawa coupling also need to be kept small enough in magnitude so as to 
not affect the stability analysis of the homogeneous case, which is possible to do with this 
setup. Such a scenario was explored for the top-quark in [51]. 

With the caveats out of the way, we now turn to the analysis. The nuclear matter in 
stars has a typical density of TVf, ss Aq CD ss (lGeV) 3 , whereas the Higgs bounce solution 
has a characteristic scale A b- From the fermion example in Sec. 3.2, we see that baryonic 
matter gives a tiny correction factor to the exponent of the decay rate, 


Se, 8 _ n ( Vu,d,A qcd\ 


(4.4) 


with y being the Yukawa coupling of a light quark commonly found in nuclear matter. 
However, even if it was the top-quark Yukawa, this exponential correction factor would 
still for all intents and purposes be one as there is such a large disparity between Aqcd 
and Ab- 

It seems fair to ask, since the correction to the exponent of the decay rate is so small, 
is the bubble formation process even sensitive to the presence of baryons. If the answer 
is no, then the probability that the Higgs vacuum decayed in the past light cone of the 
universe, P, is the same tiny number as in the homogeneous case, P = Pq. On the other 
hand, if the answer is yes, then we need to calculate the pre-exponential factor. Based on 
the analysis on Sec. 2, we find, 


P = 



A| \ 

^QCD-^O J 


g —5( A s ). 


(4.5) 


The first term in parentheses estimates the number of baryons in the visible universe. An 
equally good expression is ^PcritH^ 4 . Both quantities give the same order of magnitude for 
the number of baryons in the visible universe, Wmp. ~ 10". All cosmological parameters 
were taken from [52]. The denominator of the second term in parentheses is an estimate 
of effective surface area of a proton, with Aq^ d being the spatial radius of a proton. The 
factor of Hq in the second term appears because the proton is stable. We have assumed 
a hyper-cylindrical shape for the proton (in spacetime), as we believe it gives a better 
approximation of the pre-exponential factor. The conclusions regarding the exponent of 
the rate decay are unaffected by this assumption on the shape of a proton. By inspecting 
the following ratio, 


^(i)(A&)“ (1 ° 77)(1 °" 43)<<1 ' <4 ' 6) 

we see that the pre-exponential factor is much smaller in the heterogeneous case than in the 
homogeneous case, leading to a subdominant heterogeneous decay rate. On that note, it’s 
worth mentioning that Ref. [24] found that Hawking radiation does not classically catalyze 
Higgs vacuum decay in the SM. 
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Even though impurities from within the SM do not catalyze Higgs vacuum decay, it may 
be possible for impurities from BSM physics to do just that. The potential enhancement of 
the decay rate hinges on both the scale and the density of the new physics. In particular, 
they must both be large, relatively speaking. First of all, the stability analysis in the 
homogeneous SM assumes that there is no new physics up to the Planck scale. If the BSM 
scale is similar to or less than the instability scale in the SM, which is rather large itself 
(A j 10 11 GeV, see [3] for more information), the BSM physics may well stabilize the Higgs 
potential. In any case, a stability analysis would have to be performed. Assuming the scale 
of the BSM physics is higher than A j avoids this issue and is more in line with the spirit 
of the homogeneous SM analysis. Additionally, the results of the examples investigated in 
Sec. 3 show that the decrease of the heterogeneous action relative to the homogeneous case 
depends on the ratio, Ag SM /A^. Based on the scalar example, Eq. (3.7), GUT scale new 
physics with |A 3 1 ~ (lO 15 GeV) 2 , and the A 2 B estimated as Ai a 2 ~ (lO 1 ' GeV)“, should 

cause a 4 — 5% decrease in the heterogeneous action if their interactions with the Higgs field 
are mediated by a scalar field . 1 The decrease in action continues as Absm approaches A b, 
until Ab is reached at which point it’s not clear if the thin-wall approximation still holds 
for the scalar example of Sec. 3.1. The pre-exponential factor is governed by the density 
of impurities. In order for a modest decrease in the action to lead to an enhanced decay 
rate, as is the case for GUT scale physics, the density of impurities large enough such that 
the pre-exponential factor does not suppress the heterogeneous decay rate. 

5 Discussion and Conclusions 

It was mentioned in the introduction that in statistical mechanics systems there are typ¬ 
ically two ways to induce a phase transition, either through an impurity in the bulk or 
through a boundary of the system. Until now, we focused entirely on impurities in the 
bulk of the metastable phase. The reason is that a spatial boundary would need to be 
within our past light cone in order for it to affect the probability that the universe would 
have already decayed by now, rather than just the decay rate itself. It seems likely that 
if such a boundary was within the past light cone of the visible universe, it would have 
already been detected, which is why this scenario is not given much consideration. If there 
is a spatial boundary of the universe outside of the light cone of the visible universe, the 
boundary would affect the lifetime of the universe, but not the probability that it has al¬ 
ready decayed. Even though bubbles of true vacuum may be forming on this boundary at 
a faster rate than they would in the bulk, they would not have had enough time to reach 
our visible universe. 

4 There is an upper bound on the masses of right-handed neutrinos operating through a see-saw mech¬ 
anism of ~ IO 14 - 14 GeV from the requirement that the lifetime of the EW vacuum is longer than the age 
of the universe [53]. With this bound, the potential for right-handed neutrinos to affect the heterogeneous 
decay rate is rather small. On the other hand, the EW vacuum can be stabilized by a scalar threshold 
effect, as shown in [54]. The combination of these two effect was recently investigated in Ref. [55] with the 
conclusion that a scalar threshold can stabilize the EW vacuum even if the theory contains right-handed 
neutrinos. 
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This raises the question of what should be considered a spatial boundary. We have 
in mind a universe that is, say, topologically flat with a finite volume bounded by, say, a 
4-sphere. This is different from a universe with a non-trivial global topology. 5 Similarly, in 
extra dimensional theories, typically the extra dimensions have periodic boundary condi¬ 
tions, which would not serve as a spatial boundary. For theories without periodic boundary 
conditions there would likely a suppression of the pre-exponential factor, originating from 
making the extra dimension(s) small enough to have presently avoided detection. 

As previously mentioned, it would be interesting to relax the assumption that the 
boundary between the impurity and the metastable phase is a flat plane. This would allow 
for the study of the cases where the impurity is a string or monopole of comparable size. 
Another possibility for a future direction would be to generalize this formalism to curved 
space. 

In this work, we derived the decay rate of an unstable phase of a quantum field theory 
in the presence of an impurity in the thin-wall approximation. This derivation was based 
on the how the impurity changes the (flat spacetime) geometry relative to case of pure 
false vacuum. Two examples were given that show how to estimate some of the additional 
parameters that enter into this heterogeneous decay rate. This formalism was then applied 
to the Higgs vacuum of the SM, where baryonic matter (stars) acts as an impurity in the 
electroweak (nretastable) Higgs vacuum. We showed that the heterogeneous decay rate 
is suppressed with respect to the homogeneous case, which is to say that the conclusions 
drawn from the homogeneous case are not modified by the inclusion of baryonic matter in 
the calculation. On the other hand, we confirmed that BSM physics with a characteristic 
scale comparable to A#, can in principle lead to an enhanced decay rate contingent upon the 
change in the pre-exponential factor in the heterogeneous case relative to the homogeneous 
case, which is governed by the density of impurities in the false vacuum. 
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